Disentangling the relationship of gut microbiota, functional gastrointestinal disorders and autism: a case–control study on prepubertal Chinese boys

Emerging evidence of an altered gut microbiome in autism spectrum disorder (ASD) suggests a pathomechanism through the gut–brain axis despite the inconsistent microbiome profile reported across studies. One of the knowledge gaps in the existing ASD microbiota studies is the lack of systematic exploration of the role of comorbid functional gastrointestinal disorder (FGID) in the association of ASD and altered gut microbiome. Consequently, 92 ASD and 112 age-matched typically developing (TD) boys were profiled on general psychopathology, FGID status by Rome IV classification, and gut microbiota using 16S ribosomal RNA amplicon sequencing at the V4 hypervariable region. Compared to TD, a significant decrease in the within-sample abundance of taxa was observed in ASD, regardless of FGID status. The microbiota of ASD FGID+ and ASD FGID− clustered apart from the TD groups. The microbiota of ASD FGID+ also showed qualitative differences from that of ASD FGID− and had the highest-level Firmicutes: Bacteroidetes ratio, which was paralleled by elevated levels of anxiety and overall psychopathology. The altered gastrointestinal microbiota composition in ASD appeared to be independent of comorbid FGID. Further studies should address how FGID may mediate neuropsychiatric symptoms in ASD through inflammation along the microbiota–gut–brain axis.

Emerging evidence of an altered gut microbiome in autism spectrum disorder (ASD) suggests a pathomechanism through the gut-brain axis despite the inconsistent microbiome profile reported across studies. One of the knowledge gaps in the existing ASD microbiota studies is the lack of systematic exploration of the role of comorbid functional gastrointestinal disorder (FGID) in the association of ASD and altered gut microbiome. Consequently, 92 ASD and 112 age-matched typically developing (TD) boys were profiled on general psychopathology, FGID status by Rome IV classification, and gut microbiota using 16S ribosomal RNA amplicon sequencing at the V4 hypervariable region. Compared to TD, a significant decrease in the within-sample abundance of taxa was observed in ASD, regardless of FGID status. The microbiota of ASD FGID+ and ASD FGID− clustered apart from the TD groups. The microbiota of ASD FGID+ also showed qualitative differences from that of ASD FGID− and had the highest-level Firmicutes: Bacteroidetes ratio, which was paralleled by elevated levels of anxiety and overall psychopathology. The altered gastrointestinal microbiota composition in ASD appeared to be independent of comorbid FGID. Further studies should address how FGID may mediate neuropsychiatric symptoms in ASD through inflammation along the microbiota-gut-brain axis.
Autism spectrum disorder (ASD) is a neurodevelopmental disorder encompassing impairing neuropsychological and behavioral symptoms, as well as comorbid medical disorders 1 . Gastrointestinal (GI) symptoms are common in ASD 2 , with up to 30.5% of the clinical ASD population suffering from one functional gastrointestinal disorder (FGID) 3 , in which altered gut microbiota composition is implicated as part of the pathophysiology 4 . This is significant, as the microbiota-gut-brain axis postulated bidirectional effects between the GI and the central nervous system to explain the pathophysiology of psychiatric disorders 5 , including ASD 6 . During neurodevelopment, GI microbiota plays a crucial role during brain maturation 7 . It has been shown that within-sample diversity, that is, the alpha diversity of GI microbiota, was associated with the variability of brain connectivities in human infants to mediate emotionality 8 . Germ-free mice displayed behavioral traits that resembled autism 9 . Hence, the gut microbiome could contribute to the pathogenesis of ASD. Furthermore, as ASD children with comorbid FGID have a higher severity of rigid-compulsive behavior 10 , sensory over-responsivity 11 , and anxiety, even at prepubertal ages 12 , the altered gut microbiome in FGID may also mediate neuropsychiatric and behavioral symptoms.
The promising framework of the microbiota-gut-brain axis, coupled with technological advances in nextgeneration sequencing, has fueled the blooming of gut microbiota research in ASD over the past decade. It has been repeatedly shown that the microbiota composition in ASD children is divergent from that of typically

Results
Demographics and clinical characteristics of the study population. A total of 92 ASD and 112 TD boys were recruited and completed the study. Ages were comparable between the groups, with means of 8.43 (SD 1.54) years and 8.12 (SD 1.99) years, respectively. The maternal age was significantly older in the ASD group, while the paternal age was similar between groups. Parental educational levels were not significantly different, suggesting a similar socioeconomic status between ASD and TD 27 (Supplementary Table S1). Comparable to a local prevalence study 28 , 54.3% of the ASD participants had comorbid attention deficit hyperactivity disorder (ADHD), and 5.4% had a specific learning disorder (SLD). In addition, 43.5% of the ASD participants were taking at least one psychiatric medication, which includes methylphenidate or atomoxetine (43.5%), aripiprazole or risperidone (10.9%) and cyproheptadine (6.5%). Moreover, 62 participants were classified as FGID+ according to the Rome IV Diagnostic Questionnaires for Paediatric Functional Gastrointestinal Disorders: Parent-report Form for Children (R4PDQ), with 30 (32.6%) in the ASD group and 32 (28.6%) in the TD group. Among all the FGID+ participants, functional abdominal pain disorders (FAPD) was the most common subcategory of FGIDs (14.7%), followed by functional defecation disorders (FDD) at 12.7% and functional nausea and vomiting disorders (FNVD) at 2.9%. None of the participants was classified as having non-retentive faecal incontinence, an entity categorized under FDD other than functional constipation. Hence, participants classified as FDD were all suffering from functional constipation. The distribution of these subcategories of FGIDs is comparable between the ASD group and TD group [χ 2 = 1.619, p = 0.655], allowing a valid contrast of the gut microbiota between ASD FGID+ and TD FGID+. Body mass index (BMI) were comparable across the four groups of ASD/TD FGID+/− [p = 0.325]. Regarding psychopathology, Autism Spectrum Quotient-10 (AQ-10), Spence Children's Anxiety Scale-Parent Version (SCAS-P) and Strengths and Difficulties Questionnaire-Parent Version (SDQ) were all significantly higher in ASD, with ASD FGID+ having the highest anxiety level and overall psychopathology, followed by ASD FGID− and then the TD groups [p-values all < 0.001]. ASD FGID+ had the lowest fiber intake, while TD FGID− had the highest main nutrient. The demographic, dietary, and clinical characteristics of the four groups are summarized in Table 1. Supplementary Table S2 summarizes the comparisons of the dietary components between the four groups.
Gut microbial diversity between ASD and TD with and without FGID. A total of 12,298,104 high-quality reads and 2,645,778 sequences were obtained from the 204 participants, with an average of 12,969 sequences per sample. Moreover, 3162 features were identified in all samples at a 99% similarity level, with an average of 836 features observed in each sample. Alpha diversity analysis revealed that the within-subject species diversity and abundance were significantly lower in ASD, as measured by both Faith's Phylogenetic Diversity (Faith's PD) [p adjusted = 0.001] and the Chao1 Index [p adjusted < 0.001] ( Table 2 and Supplementary Fig. S1). For beta diversity, the overall gut microbial composition of ASD and TD was significantly different, as measured by the Bray-Curtis dissimilarity [p adjusted = 0.001], unweighted UniFrac distance metrics [p adjusted = 0.001], and weighted UniFrac distance metrics [p adjusted = 0.001] ( Table 3 and Supplementary Fig. S2).
Stratifying both ASD and TD into groups of FGID+ or FGID− for pairwise comparison showed that regardless of FGID status, the ASD population had a decrease in the abundance of taxa, measured by the Chao1 index, compared to the corresponding TD groups [FGID−: p adjusted < 0.001; FGID+: p adjusted 0.04], whereas phylogenetic diversity, as measured by Faith's PD, was lower in ASD in the FGID− [p adjusted < 0.001] but comparable between ASD and TD with FGID. When contrasting ASD with and without FGID, although the Chao1 index was comparable, ASD FGID+ participants had higher Faith's PD [p adjusted 0.032]. Conversely, the alpha diversity was similar between TD with and without FGID ( Table 2    www.nature.com/scientificreports/

Differences in taxa between ASD and TD with and without FGID.
Further analyses compared the relative abundance of taxa between ASD and TD. Firmicutes, Bacteroidota, Actinobacteriota, Proteobacteria and Desulfobacterota constituted the most abundant phylum among the participants. Significant variations were observed, and comparisons of the relative abundance at the phylum level between the groups are summarized comprehensively in Supplementary Table S7 and Supplementary Fig. S4. Notably, the Firmicutes: Bacteroidetes (F:B) ratio was observed to be higher in ASD than TD [M ASD = 10.9, M TD = 5.78, p = 0.003]. Furthermore, ASD FGID+ had the highest F:B ratio, following a descending trend in ASD FGID− and TD FGID+ to TD FGID− ( Table 4, Fig. 2). The genus-level analysis identified 19 genera that were significantly different between ASD and TD (Supplementary Table S8 and Fig. 3A,F). Discriminative genera were identified in each of the four participant groups (Fig. 3B,G). Compared to TD FGID-, ASD FGID-showed a significant increase in Bifidobacterium, Dorea, and Blautia and a decrease in Collinsella, Bacteroides, Alistipes, Parabacteroides, Clostridium sensu stricto 1, Sutterella, and uncultured Eggerthellaceae. ASD FGID+ was distinctive from ASD FGID-and TD FGID+, with a significant increase in Ruminococaceae (family), Dialister, Fusicatenibacter, Turicibacter, and a decrease in Phascolarctobacterium (Table 5, Fig. 3C-E). The comparison between TD FGID+ and TD FGID− is summarized in Supplementary Table S9.
As the ASD groups were characterized by a lower alpha diversity and a higher F:B ratio, correlations between these indices with clinical parameters of AQ-10, SCAS-P and SDQ were also explored. It was observed that the lower alpha diversity and higher F:B ratio were significantly associated with higher AQ-10 and SDQ scores (Supplementary Table S10).

Discussion
In this case-controlled sample of prepubertal ASD boys and age-and gender-matched TD, the altered microbiota composition was observed in ASD over TD, considering the confounding effect of comorbid FGID. The gut microbiome of ASD comorbid with FGID differed from ASD without FGID by phylogenetically distinct taxa that were of low relative abundance. Furthermore, autistic symptoms (AQ-10), anxiety (SCAS-P) and overall psychopathology (SDQ) levels were highest among the ASD FGID+ and explorative analyses also showed that lower alpha diversity and microbiota composition at the phylum level (F:B ratio) correlated with higher SDQ and AQ-10. Our results thus corroborated the reported potential correlations between microbiota, GI and neuropsychiatric symptoms in ASD [10][11][12] . Notwithstanding the rigorous attempts made to control for possible confounders and optimize the characterization of psychopathology and FGIDs in this study, our findings should be interpreted from the wider perspective of mixed results in published studies. It is also worth noting that, a higher or lower alpha diversity do not have a direct implication on the host's health status and should be interpreted as a consequence of the ecological mechanism that drives the observations 29 . Therefore, the compositional changes of microbiota in ASD in the present study, regardless of the directionality, should prompt further investigation on the underlying biological cause. Further, the degree and pattern of disease-associated microbiota changes vary across different entities, with enrichment of a small quantity of pathogens, depletion of health-related taxa and gross compositional changes being reported 30 , and little is known in this area for ASD. www.nature.com/scientificreports/ Alteration of gut microbiota composition in ASD: merely a confounding effect of FGID or an independent signature? There were inconsistent results across alpha diversity analyses in ASD, where a decrease, increase, and most commonly, similar richness and diversity were reported 31 . Our results showed that, when excluding FGID, a homogenous decrease in alpha diversity across the two indices measured was observed in ASD. Also, beta diversity showed the microbiota in ASD are distinct from TD in subjects without FGID across the three indices measured. Collectively, these findings suggest that the altered abundance, diversity, and structure of microbiota in ASD is not merely a result of the higher prevalence of comorbid FGID. Within our sample, several genera were found to be significantly different and discriminative between ASD and TD without FGID. However, some of the results conflicted with a recent systematic review 13 . For example, Bifidobacterium [32][33][34][35] and Blautia 36,37 were both reported to have a lower relative abundance in ASD, while an increase in Collinsella 38 was observed. Mixed results (both increase and decrease) were also reported for Dorea, Bacteroides, Alistipes, Oarabacteroides, and Sutterella 13 . Here, several possible reasons are discussed for the inconsistencies in methodological heterogeneity among the studies included in the systematic review. First, some studies quoted either did not separate ASD participants with and without GI symptoms or FGID prior to analysis 33,38 , or included only ASD participants with comorbid GI symptoms 32,36,37 . Second, the gender distribution between ASD and the controls was unequal for some studies [33][34][35][36] , with females being over-represented in the control group. This may confound the results, given that Bifidobacterium were reported to be higher in females 39 . Third, compared to the present study, the sample sizes of these previous studies were all relatively smaller, in the range of six to 40 ASD cases. Despite so, the result of increase in abundance of Bifidobacterium in ASD FGID− relative to TD FGID− warrants further investigation, as the genus is frequently cited as beneficial 40 .
Altered gastrointestinal microbiota composition is implied in the pathophysiology of ASD. It has been recently shown that host genetics and gut microbiota interact to regulate immune and metabolic pathways relevant to ASD 41 , and various pathways have been postulated on how the microbiome can exert its influence on the host's brain. These include the regulation of gut permeability through tight junction protein, the production of shortchain fatty acids (SCFAs), and the metabolism of tryptophan and inflammatory cytokines 7 . Within the genera that discriminated between FGID-ASD and TD from the present study, Alistipes, Bifidobacterium, Blautia, and Collinsella are all fermenters that produce SCFAs, while Dorea and Sutterella could regulate the mucosal metabolism to affect gut permeability 42,43 . With a lower relative abundance of Bacteroides in the ASD FGID− group observed, it has been shown that oral treatment with Bacteroides fragilis corrected ASD-related behavior in mice by altering gut permeability and metabolic pathways 44 .

Is ASD FGID+ a distinctive subgroup in the gut microbiome profile?
Within the smaller sample of FGID+ participants, the lower Chao1 index and beta diversity analysis indicated that the microbiota of ASD FGID+ had a lower taxa abundance and were clearly clustered apart from that of TD FGID+. Although there are distinct pathophysiologies of each FGID subtype 45 , given the comparable distribution of FGID subtypes in the ASD and TD groups in this study, our observations are unlikely to be confounded by the biological factors unique to individual subtypes of FGID. Hence, the altered gut microbiota composition in ASD FGID+ did not appear to be solely attributable to FGID, but possibly reflect an additive effect of both ASD-associated and FGIDassociated changes in the microbiome. Furthermore, with higher Faith's PD in the ASD FGID+ group relative to ASD FGID−, and the two groups clustered only in the measurement of unweighted UniFrac distance but not weighted UniFrac and Bray Curtis dissimilarity, it appears the microbiota of ASD FGID+ was phylogenetically more diverse, and these phylogenetically distinct taxa that presented in low relative abundance were setting ASD FGID+ apart from ASD FGID−. Comparatively, the distinction between ASD FGID+ and ASD FGID− was not as pronounced as when ASD was compared with TD. Hence, the microbiota of ASD with and without FGID could lie along a continuum within the altered gut microbiota composition of ASD.
As an increase in the F:B ratio was suggested to be linked with a proinflammatory state, there was a gradient of increase of the F:B ratio from TD FGID−, ASD FGID−, and TD FGID+ to the highest in the ASD FGID+ participants. Intriguingly, this gradient ratio also paralleled the clinical observation that the ASD FGID+ group had the highest levels of anxiety and overall psychopathology. Thus, the correlation between FGID and neuropsychiatric symptoms in ASD could be mediated by a proinflammatory state along the microbiota-gut-brain axis that warrants further exploration.
Among the genera that distinguished ASD FGID+ from both ASD FGID− and TD FGID+, Fusicatenibacter, Dialister, and Phascolarctobacterium are producers of SCFAs. Our results also showed that ASD FGID+ was enriched with Turicibacter, the abundance of which was found to have a unique reciprocal positive effect on the gut luminal level of serotonin (5-HT) 46 . Turicibacter, together with other spore-forming bacteria dominated by the clostridia class, promote the gut production of 5-HT by inducing the rate-limiting enzyme tryptophan hydroxylase-1, while 5-HT conversely promotes the competitive colonization of Turicibacter in the intestine. Gutderived 5-HT, which constitutes > 90% of the human body's 5-HT 47 , is an important modulator of gastrointestinal motility 48 , and an altered level of 5-HT in the gastrointestinal system was implicated in the pathophysiology of FGIDs 49 . Thus, the increase in the relative abundance of Turicibacter observed in our ASD FGID+ participants may illuminate the pathophysiology of GI symptoms in ASD individuals. Furthermore, the 5-HT produced in the gut is stored in circulating platelets and may result in hyperserotonaemia, a phenomenon observed in 30% of ASD cases 50 . While the phenomenon remained enigmatic in its clinical correlates and causes, our results revealed the altered gut microbiota in ASD may provide clues to its occurrence.

Future directions and limitations
Although the present study is insufficient to infer or elucidate the underlying pathophysiological mechanism of microbiota changes in ASD, the clarification of their relationship with FGID serves as a starting point to further understand the causes and effects in young children. Direct measurement of the functional output of the gut microbiome through a multi-omics approach, including metatranscriptomics, metabolomics, and metaproteomics, could further reveal the downstream effect of microbiota to infer the possible pathomechanism 51 . In our present study, a functional analysis such as PICRUSt 52 was not conducted for two main reasons. First, the main objective was to verify the change in the microbiota composition in ASD independent of comorbid FGID rather than to propose any possible pathomechanism. Second, there is an inherent limitation of 16 s amplicon sequencing in terms of the taxonomic assignment beyond the genus level and thus the accuracy of functional analysis 53 . Beyond that, the cross-talk between the microbiome and the host's genetic expression is also worth investigating as a potential pathogenesis pathway 54 . ASD is heterogeneous in its severity, which is multi-dimensionally captured by the core ASD symptoms, comorbid psychopathology and functional level 55 , and whether gut microbiota has a role in the severity of core symptoms and comorbid psychopathology in ASD is scarcely addressed. Our explorative analyses showed that microbiota composition correlates with these clinical aspects of ASD, and within the ASD FGID− subgroup, microbiota was distinctive between those with and without ADHD. A systematic investigation by profiling the ASD symptom domains and comorbid psychopathology comprehensively, and comparisons with non-ASD individuals with similar psychopathology, would shed light on the potential mediating and moderating role of gut microbiota in the diverse clinical phenotypes of ASD. The medications' effect on gut microbiota composition in ASD is also to be studied by considering individual medication and the dosage effect.
Several methodological limitations should be considered when interpreting the findings of this study. First, while a stringent sampling frame was implemented to minimize the confounding effects of puberty and gender, this would also limit the generalizability of the results. Second, some stool samples stored overnight at the participants' domestic refrigerator at 4 °C were not ideal and could introduce spurious results. Third, the dietary assessment captured merely an overall food composition, which lacked information on specific food types and allergy, dietary preference, and history of dietary intervention. Comprehensive profiling of these would help uncover potential mediating and moderating effects on the gut microbiota by diet, as suggested in a recent study 56 . Fourth, the sample size of FGID+ was relatively small, which is a heterogeneous group of disorders with divergent pathophysiology. The different stool transit and fermentation time in constipation and diarrhoea may result in different bacterial composition. Hence, categorizing different FGIDs as a single group in the analyses has inevitably precluded the study of biodiversity across the different types of FGIDs. A sufficiently powered study with the full range of FGIDs would be needed in the future to investigate their specific associations with the gut microbiome in ASD. Also, there was a lack of dimensional measurement of GI symptoms to consider the spectrum of varying FGID severity and overlapping symptoms 57 . Finally, the evaluation of the anxiety level, psychopathological severity, and FGID status was respondent-based, which may introduce response biases.

Conclusion
The presence of altered gut microbiota composition was associated with ASD after accounting for the presence of FGID. A small quantity of phylogenetically distinct taxa characterized the differences in the gut microbiome between ASD with and without FGID. Further understanding of the pathomechanism in the microbiota-gut-brain axis in ASD would benefit from functional analyses of the gut microbiome and the utilization of a multi-omics approach in combination with network analysis.

Recruitment and clinical assessment of research participants.
In this case-control study, Chinese ASD and TD boys aged between 4 and 11 years old were recruited. Other than the presence of FGID in the FGID+ groups, all participants did not have other major illnesses such as organic GI or neurological disorders. The age range of the participants was confined to 4-11 years old to minimize the effect of puberty 58 , while ensuring the structural stability of the microbial community, which is generally reached at 3 years of age 59 . Only boys were studied, as the effect of microbiota on the brain could be gender-dependent 8  Community-based TD children with reported good general physical and mental health were recruited from public schools. Participants and their parents were informed of the details of the study, and written informed consent was obtained from the parents according to the Declaration of Helsinki. Clinical records were reviewed for ASD boys, and demographic data were obtained from the parents of all the participants. Body height and weight were measured to calculate their BMI. All parents were asked to complete the following questionnaires: (1) The R4PDQ, an instrument designed to assess FGID in children in both clinical and research settings according to the Rome IV diagnostic criteria 45 . From the three broad categories of FGID, namely FNVD, FAPD and FDD, 10 specific disorders could be identified. The R4PDQ was validated against a gastroenterologist's diagnosis and demonstrated diagnostic accuracy as a respondent-based classification instrument 60 64 . Parents completed the questionnaire on children's eating behaviors and dietary habits to obtain information on their main nutrient intake (the amount of milk, fruit, vegetables, cereal, meat, fish, beans and egg consumed), fiber intake (the amount of fruit and vegetables consumed), and unhealthy intake (the amount soft drink and fried food consumed, habit of eating at fast food and fatty meat) over the past seven days. Details of the dietary questionnaire and the algorithm for calculating the three composite scores is listed in the Supplementary Information. The exclusion criteria were: (1) A score above the cut-off (≥ 6) on the AQ-10 for the TD group 61 ; (2) Intellectual disability and known history of organic GI disorder based on the medical record; (3) The use of antibiotics within a month before stool sample collection; and (4) The regular daily use of probiotics within a month prior to stool sample collection. Bioinformatics and data processing. Data processing was performed using the QIIME 2 bioinformatics software package 65 . Raw paired FASTQ sequences files were demultiplexed using the q2-demux plugin and controlled for sequence quality with the processing pipeline DADA2 (q2-dada2). The processing of reads with the dada2-paired method was preceded by trimming the first five bases and truncating at 300 bases in all forward and reverse reads. Reads were filtered, dereplicated and merged to generate the following: (1) A feature table containing counts of 16S rRNA gene sequence variants for each sample; (2) A representative sequence table containing the sequences for each unique feature. A taxonomy was then assigned to each feature sequence against the Silva-138 database using the classify-consensus-vsearch implemented in the q2-feature-classifier plugin. Feature tables were further filtered by removing sequences identified as archaea, mitochondria, or chloroplasts. Following the taxonomy assignment, the feature sequences were aligned with the q2-alignment plugin MAFFT. An unrooted tree using the FastTree method and a rooted tree using the midpoint rooting method were generated with the q2-phylogeny plugin. The QIIME 2 output files, including a feature table, taxonomy, and unrooted tree, were imported to R software (Version 4.0.3) to perform further data processing and analysis. Feature tables were rarefied to 6808 sequences per sample.
Statistical analysis. Demographics of the participants were compared using the Chi-square test for categorical data, the Student t-test, or the Wilcoxon rank-sum test for continuous data, depending on the normality. Both ASD and TD participants were stratified into FGID+ and FGID− subgroups according to the Rome IV diagnostic criteria. An overall comparison between ASD and TD was conducted before subgroup analyses pertaining to our specific study objectives. First, we compared ASD FGID− vs. TD FGID− and ASD FGID+ vs. TD FGID+ to test if ASD has altered microbiomes both in the presence and absence of FGID. Second, we compared ASD FGID+ vs. ASD FGID− for disparities in the microbiome between ASD with and without FGID. TD FGID+ vs. TD FGID− was also conducted as a reference. Several indices were computed using the phyloseq R package to compare the compositionality of the gut microbiome between groups, with each index having its own approach to measure the ecological structure. First, to contrast the compositional complexity between ASD and TD (i.e. within-sample diversity or alpha diversity), the Chao1 index 66 and Faith's PD 67 were estimated and compared between groups with adjusted multivariate linear regression analyses. The Chao1 index is an abundance-based estimator of species richness, while Faith's PD considers phylogenetic distance when estimating the diversity. Subsequently, to test for the taxonomical differences between ASD and TD (i.e. between-sample diversity or beta diversity), Bray-Curtis dissimilarity 68 , unweighted, and weighted UniFrac 69 were used. Bray-Curtis dissimilarity quantifies dissimilarity based on abundances of shared taxa, while UniFrac takes the phylogenetic tree into consideration for compositional dissimilarity, with unweighted UniFrac calculated on the mere presence or absence of the taxa and weighted UniFrac calculated on the disparity in the proportional abundance of the taxa. Hence, unweighted UniFrac is more sensitive to differences in low-abundance features. The permutational Multivariate Analysis of Variance (PERMANOVA) was performed on these three indices using the adonis function imple- www.nature.com/scientificreports/ mented in the vegan R package with 1000 permutations, which were then visualized using principal coordinates analysis (PCoA). All testing on alpha diversity and beta diversity estimates were performed with adjustments to the usage of psychiatric medications and dietary information as covariates.
After the compositional structures, variations in the relative abundance between the ASD and TD groups at the phylum and genus levels were compared using the Wilcoxon rank-sum test, assuming the data of relative abundance were not normally distributed. The F:B ratio, which is composed of the two most abundant phyla of the human gut microbiome, was measured and compared. An increased F:B ratio was reported in ASD 38 , and inflammatory bowel disease 70 and was proposed to be related to inflammatory status 71 . P-values were corrected for the false discovery rate (FDR) upon multiple comparisons using the Benjamini-Hochberg method. Finally, the Linear discriminant analysis Effect Size (LEfSe) analysis was performed using the online Galaxy Hutlab LEfSe platform 72 . Kruskal-Wallis tests within the LEfSe were implemented to analyze all features to identify the differentially abundant features that distinguish the groups, with the threshold of statistical significance at 0.05 and the logarithmic linear discriminant analysis (LDA) score for discriminative features at 2.0.

Data availability
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.